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ABSTRACT 

A toy model for gamma-ray burst supernovae (GRB-SNe) is discussed by considering 
the effects of screw instability of magnetic field in black hole (BH) magnetosphere. The 
screw instability in the Blandford-Znajek (BZ) process (henceforth SIBZ) can coexist with 
the screw instability in the magnetic coupling (MC) process (henceforth SIMC). It turns 
out that both SIBZ and SIMC occur inevitably, provided that the following parameters are 
greater than some critical values, i.e., (i) the BH spin, (ii) the power-law index describing 
the magnetic field at the disk, and (iii) the vertical height of the astrophysical load above 
y—i ■ the equatorial plane of the rotating BH. The features of several GRBs are well fitted. In 
^ our model the durations of the long GRBs depend on the evolve time of the half-opening 

angle. A small fraction of energy is extracted from the BH via the BZ process to power a 
GRB, while a large fraction of energy is extracted from the BH via the MC process to power 
an associated supernova. In addition, the variability time scales of tens of msec in the light 
\Q ', curves of the GRBs are fitted by two successive flares due to SIBZ. 
O ; 

Subject headings: accretion, accretion disks - black hole physics - gamma rays: bursts- 
supernovae: general — magnetic fields — instability 
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1. INTRODUCTION 

> ■ 

• i-H , 

Recently, observations and theoretical considerations have linked long-duration GRBs with ultra- 
■ bright Type Ib/c supernovae (SNe; Galama et al. 1998, 2000; Bloom et al. 1999). The first candidate 
was provided by SN 1998bw and GRB 980425, and the recent HETE-II burst GRB 030329 has greatly 
enhanced the confidence in this association (Stanek et al. 2003; Hjorth et al. 2003). Extremely high 
energy released in very short time scale suggests that GRBs involve the formation of a black hole (BH) 
via a catastrophic stellar collapse event or possibly a neutron star merger, implying that an inner engine 
could be built on an accreting BH. 

Among a variety of mechanisms of powering GRBs the Blandford-Znajek (BZ) process (Blandford 
& Znajek, 1977) has its unique advantage in providing "clean" (free of baryonic contamination) energy 
by extracting rotating energy from a BH and transferring it in the form of Poynting flow in the outgoing 
energy flux (Lee et al. 2000, hereafter LeeOO; Li 2000c, hereafter LiOO). 

Not long ago Brown et al. (2000, hereafter BOO) worked out a specific scenario for a GRB-SN 
connection. They argued that the GRB is powered by the BZ process, and the SN is powered by the 
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magnetic coupling (MC) process, which is regarded as one of the variants of the BZ process (Blandford 
1999; van Putten 1999; Li 2000b, 2002; Wang, Xiao & Lei 2002, hereafter W02). It is shown in BOO 
that about 10 53 ergs are available to power both a GRB and a SN. However, they failed to distinguish 
the fractions of the energy for these two objects. 

More recently, van Putten and his collaborators (van Putten 2001; van Putten & Levinson 2003, 
hereafter P03) worked out a poloidal topology for the open and closed magnetic field lines, in which 
the separatrix on the horizon is defined by a finite half-opening angle. The duration of a GRB is set 
by the lifetime of the rapid spin of the BH. It is found that GRBs and SNe are powered by a small 
fraction of the BH spin energy. This result is consistent with observations, i.e., duration of GRBs of 
tens of seconds, true GRB energies distributed around 5 x 10 50 ergs (Frail et al. 2001), and aspherical 
SNe kinetic energies of 2 x 10 51 ergs (Hoflich et al. 1999). 

Very recently, Lei et al. (2005, hereafter Lei05) proposed a scenario for GRBs in type Ib/c SNe 
invoking the coexistence of the BZ and MC processes. In Lei05 the GRB is powered by the BZ process, 
and the associated SN is powered by the MC process. The overall time scale of the GRB is fitted by 
the duration of the open magnetic flux on the horizon. 

Besides the features of high energy released in very short durations most GRBs are highly variable, 
showing very rapid variations in flux on a time scale much shorter than the overall duration of the 
burst. Variability on a time scale of milliseconds has been observed in some long bursts (Norris et al. 
1996; McBreen et al, 2001; Nakar and Piran, 2002). Unfortunately, the origin of the variations in the 
fluxes of GRBs remains unclear. In this paper we intend to discuss the mechanism for producing the 
variations in the fluxes of GRBs by virtue of the screw instability in BH magnetosphere. 

It is well known that the magnetic field configurations with both poloidal and toroidal components 
can be screw-unstable. According to the Kruskal-Shafranov criterion the screw instability will occur, if 
the toroidal magnetic field becomes so strong that the magnetic field line turns around itself once or 
more (Kadomtsev 1966; Bateman 1978). 

Some authors have addressed the screw instabilities in BH magnetosphere. Gruzinov (1999) argued 
that the magnetic field with a bunch of closed field lines connecting a Kerr BH with a disk can be screw- 
unstable, resulting in the release of magnetic energy with the flares at the disk. Li (2000a) discussed 
the screw instability of the magnetic field in the BZ process, leading to a stringent upper bound to the 
BZ power. Wang et al. (2004, hereafter W04) studied the screw instability in the MC process. They 
concluded that this instability could occur at some place away from the inner edge of the disk, provided 
that the BH spin a* and the power-law index n for the variation of the magnetic field on the disk are 
greater than some critical values. 

In this paper we attempt to combine the screw instability of the magnetic field with the coexistence 
of the BZ and MC processes. To facilitate the description henceforth we refer to the screw instability of 
the magnetic field occurring in the BZ and MC processes as SIBZ and SIMC, respectively. It is shown 
that both SIBZ and SIMC can occur, provided that the following parameters are greater than some 
critical values: (1) the BH spin, (2) the power-law index describing the variation of the magnetic field 
at the disk, and (3) the vertical height of the astrophysical load above the equatorial plane of the Kerr 
BH. 

The features of several GRB-SNe are well fitted in our model. (1) The overall duration of the GRBs 
is fitted by the evolution of the half-opening angles. (2) The true energies of several GRBs are fitted 



- 3- 



by the energy extracted in the BZ process, and the energies of associated SNe are fitted by the energy 
transferred in the MC process. (3) The variability time scales of tens of msec in the light curves of 
several GRBs are fitted by two successive flares due to SIBZ. 

This paper is organized as follows. In § 2 we derived a criterion of SIBZ based on the Kruskal- 
Shafranov criterion and some simplified assumptions on the remote load. In § 3 we discuss the time 
scale and energy extraction from a Kerr BH in the context of the suspended accretion state. In § 4 we 
propose a scenario for the origin of the variation in the light curves of GRBs based on the flares arising 
from SIBZ. Finally, in § 5, we summarize the main results and discuss some issues related to our model. 
Throughout this paper the geometric units G — c — 1 are used. 



2. SCREW INSTABILITY IN BH MAGNETOSPHERE 

In W04 the criterion of SIMC is derived based on the following points: (1) the Kruskal-Shafranov 
criterion, (2) the mapping relation between the angular coordinate on the BH horizon and the radial 
coordinate on the disk, and (3) the calculations of the poloidal and toroidal components of the magnetic 
field at the disk. The criterion of SIBZ can be derived in an analogous way. However, the BZ process 
involves unknown astrophysical loads, to which both the mapping relation and the calculations for the 
poloidal and toroidal components of the magnetic field are related. In order to work out an analytical 
model we present some simplified assumptions as follows. 

(1) The magnetosphere anchored in a Kerr BH and its surrounding disk is described in Boyer- 
Lindquist coordinates, in which the following Kerr metric parameters are involved (MacDonald and 
Thorne 1982, hereafter MT82). 

£ 2 = (r 2 + a 2 ) 2 - a 2 A sin 2 6, p 2 = r 2 + a 2 cos 2 9, 

A = r 2 + a 2 -2Mr, w= (E/p)sin0, (1) 
a = pVA~/z. 

(2) The remote load is axisymmetric, being located evenly in a plane with some height above the 
disk. In Figure 1 the open magnetic field lines connect the BH horizon with the load. The symbol Lbz 
and H c represent the critical field line and the height of the remote load above the equatorial plane for 
the occurrence of SIBZ, respectively. 

(3) In Figure 1 the radius r s is the critical radius of SIMC, which is determined by the criterion of 
the screw instability given in W04, 

{2itm D /L MC )B p D /B T D <l. (2) 

In equation (2) Lmc is the critical length of the poloidal field line for SIMC, and B P D and Bp are 
the poloidal and toroidal components of the magnetic field on the disk, respectively, and w D is the 
cylindrical radius on the disk and it reads 

w D = H D /p D = {MxLVi + aU~ 2 X ~i + 2a 2 £- 3 X - 6 s . (3) 
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where Xms = \fr ms /M is defined by Novikov & Thorne (1973) in terms of the radius of innermost stable 
circular orbit (IS CO). 

(4) The angle 0$ in Figure 1 is the half-opening angle of the magnetic flux tube on the horizon, 
which is related by the mapping relation between the angular coordinate on the BH horizon and the 
radial coordinate on the disk as follows (Wang et al., 2003, hereafter W03), 

cos9 — cos9 L = G (a*;£,n) d£, (4) 

where 



* ~ 2y/(l + alx~i + 2a*x£) (1 - + «1 ' 

In equations (4) and (5) £ = r/r ms is defined as a radial parameter in terms of r ms , and n is a power-law 
index for the variation of the poloidal magnetic field at the disk, i.e., 

B p D oc r n - (6) 

(5) The suspended accretion state is assumed due to the transfer of angular momentum from the 
BH to the disk (van Putten & Ostriker, 2001). 

Analogous to equation (2) the criterion for SIBZ can be expressed as 

(2nR/L BZ )B p L /Bl<l, (7) 

where Lbz is the critical length of the poloidal field line for SIBZ, and B P L and B\ are the poloidal 
and toroidal components of the magnetic field on the remote load, respectively, and R is the cylindrical 
radius of the remote load with respect to the symmetric axis of the BH. 

The toroidal field component B\ can be expressed by Ampere's law, 

B T L = 2I L /R, (8) 
where II is the electric current flowing in the loop KMM'K' in Figure 1 and it reads 

h = VPbz/Z l . (9) 

The quantities Pbz and Zl in equation (9) are the BZ power and the load resistance, respectively. 
The BZ power has been derived in W02 as follows, 

= * f' \ {l a k) tZ (10) 

Jo 2 - (1 - q) sin 9 
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where q = ^/l — a 2 is a parameter depending on the BH spin, and k = VLp/£l H is the ratio of the 
angular velocity of the field lines to that of the BH horizon. The quantity P is defined by 

P = (B P H ) 2 M 2 « 6.59 x 10 50 x Bf 5 (M/M ) 2 erg ■ s _1 , (11) 

where B i5 represents the magnetic field at the BH horizon in terms of 10 15 gauss . 

MT82 argued in a speculative way that the ratio k will be regulated to about 0.5 by the BZ process 
itself, which corresponds to the optimal BZ power with the impedance matching. Taking the impedance 
matching into account, we have the remote load resistance equal to the horizon resistance, and they 
read 

AZ L = AZ H = R H p^-, (12) 

where Rh = 47T = 377ohm is the surface resistivity of the BH horizon (MT82). Thus we have Z L and 
Z H expressed as 

Zl = Z H =f- Ra £& (13) 

Jo 2,71ZU H J Q w H 

Incorporating equations (8) — (13), we can calculate B\ in terms of the cylindrical radius R. On 
the other hand, the poloidal magnetic field B^ at the radius R of the load can be determined by the 
conservation of the magnetic flux, i.e., 



B p H 2ixw H p H de = B[2irRdR. (14) 



From equation (14) we have 



D P _ B H W hPh d0 

~ R~dR m { } 

Assuming that the height of the planar load above the equatorial plane of the Kerr BH isif , we 
have an approximate relation between the angle 9 and the radius R as follows, 

t&n9 = R/H. (16) 
Substituting equation (16) into equation (15), we have 

P B p h w hPh cos 2 9 
Bl = HR • (17) 

Incorporating equations (8) and (17) with the criterion (7), we have 
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7iB^w H p H sin 6 S cos 2 6 S < l _ 
H^P BZ /Z L ' 

where the relation sin 6s = R/Lbz is used. The criterion (18) implies that SIBZ will occur, provided 
that the height of the load is greater than the critical height, i.e., H > H c , and we have 

h ^ H /M = *Bfa H p sinW ^ 

As argued in W04 the angle 6s can be determined by the criterion of SIMC with the mapping 
relation between the BH horizon and the disk, and it is a function of the BH spin a* and the power-law 
index n, i.e., 6s (a*, n). Inspecting equations (1), (10) and (19), we find that h c is a dimensionless 
parameter also depending on the parameters a* and n, i.e., h c = h c (a*,n). By using equations (2) and 
(19) we have the contours of 6s (a*, n) and h c (a*, n) in a* — n parameter space as shown in Figure 2. 

Inspecting Figure 2, we find the following features of the contours: 

(1) The values of 6 S increases and those of h c decreases with the increasing n for the given BH spin 
a*, respectively. 

(2) The values of 6 S increases and those of h c decreases with the increasing a* for the given n, 
respectively. 

(3) Both SIMC and SIBZ will occur, provided that the parameters a*, n and h c are greater than 
some critical values. As shown in Figure 2c the shaded region indicates the value ranges of a* and n in 
which both 6s > and h c > 100 are constrained. Thus the occurrence of SIMC and SIBZ is guaranteed 
by the value ranges of a* and n in the shaded region. 



3. TIME SCALE OF A GRB AND ENERGY EXTRACTION FROM A ROTATING 

BH 

There are several scenarios to invoke the BZ process for powering GRBs, and the main differences 
among these scenarios lie in the environment of a spinning BH and the approaches to the duration of a 
GRB. These scenarios are outlined as follows. 

Model I: In LeeOO the energy is extracted magnetically from a rotating BH without disk, and the 
duration of a GRB is estimated as the time for extracting all rotational energy of the central BH via 
the BZ process. 

Model II: It is argued that the energy is extracted magnetically from a rotating BH with a transient 
disk, and the duration of a GRB is estimated as the time for the disk plunged into the BH (Lee & Kim 
2002; Wang et al. 2002). 

Model III: In LiOO the energy is extracted magnetically from a rotating BH with a stationary 
torus in the state of suspended accretion, and the duration of a GRB is estimated roughly as the time 
for extracting all rotational energy of the central BH via the BZ process. 
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Model IV: In BOO the energy is extracted magnetically from a rotating BH with a non-stationary 
disk in the state of suspended accretion, and the duration of a GRB is determined by the presence of 
the disk. 

Model V: In P03 the energy is extracted magnetically from a rotating BH with a torus in the 
state of suspended accretion, and the duration of a GRB is determined by the instability of the disk. 

Model VI: In Lei05 the energy is extracted magnetically from a rotating BH with a thin disk 
in the state of suspended accretion, and the duration of a GRB is determined by the lifetime of the 
half-opening angle. 

This Model: It is a modified version of Model VI, in which the effects of SIMC and SIBZ are 
taken into account. Compared with Model VI some features and advantages are given as follows. 

(1) The magnetic field configuration in Model VI is built based on the conservation of the closed 
magnetic flux connecting the BH with the disk with the precedence over the open magnetic flux, resulting 
in the closed field lines connecting the half-open angle 9bz at horizon with the disk extending to infinity. 
In this Model, however, the closed field lines are confined by SIMC within a region of a limited radius 
r s (about a few Schwarzschild radii as shown in Table 4), which is consistent with the collapsar model 
for GRBs-SNe. 

(2) As argued in the next section, the variability time scales of the light curves of GRBs are 
modulated by two successive flares due to SIBZ. 

The main features of the above models for GRBs are summarized in Table 1. 

In Model VI the duration of a GRB is regarded as the lifetime of the half-opening angle 9 BZ , which 
is based on the evolution of the rotating BH. The same procedure can be applied to this model except 
that the angle 6bz is replaced by 9s arising from SIMC. It is found that the characteristics of the BH 
evolution in this model are almost the same as given in Model VI as shown by a* — n parameter spaces 
in Figure 3. 

For several GRB-SNe, the observed energy E 7 and duration T 90 can be fitted by adjusting the 
parameters n and E>i 5 . The energy Esn can be predicted as shown in Table 2, where the values of n, 
B 15 and E S n fitted to five GRBs invoking Model VI and this model are shown in the left and right 
sub-columns, respectively. 

The fractions of extracting energy from a rotating BH via the BZ and MC processes are defined 
respectively as Jbz and fuci and they read 

E BZ E M c (9n , 

t^BZ + i^MC &BZ + tliMC 

where Ebz and Emc are the energies extracting in the BZ and MC processes, respectively. 

rtBz rtBz 
Ebz = / PBzdt, E MC = / Pucdt. (21) 
Jo Jo 



In equation (21) t B z is defined as the lifetime of the angle 9s, which can be calculated by the same 
procedure given in Lei05. The MC power in equation (21) is expressed as (W04) 
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Pmc/Po = 2al 



J 

J 8. 



r L 13(1-13) sin 3 9d6 
6s 2-(l-g)sin 2 fl' 



(22) 



where the parameter (3 = VLp/Vt H = Q D /Q H is the ratio of the angular velocity of the magnetic field 
lines to that of the BH. 

In Model VI the BH spin af RB corresponds to the time when the half-opening angle 9 B z = 0, and 
the BZ power P BZ = 0, while af N corresponds to the MC power P MC = 0. In this model af RB and af N 
have the same meanings as given in Model VI except that 6 BZ is replaced by 6 S - 

As shown in Table 3, the values of af RB , f BZ and f MC fitted to five GRBs invoking Model VI and 
this model are shown in the left and right sub-columns, respectively. 

Inspecting the data in the left and right sub-columns in Table 2 and Table 3, we find that f BZ 
and Jmc i n this model are greater and less than the counterparts in Model VI, respectively. This result 
arises from the effects of SIMC and SIBZ: the half-opening angle 6$ in this model is greater than the 
half-opening angle 6 B z in Model VI as argued in W04. 



As is well known, the bursts are divided into long and short bursts according to their Tgo. Most 
GRBs are highly variable, showing 100% variations in flux on a time scale much shorter than the overall 
duration of the burst. The bursts seem to be composed of individual pulses, with a pulse being the 
"building block" of the overall light curve. The variability time scale <5tis much shorter than the GRBs' 
duration T go , the former is more than a factor of 10 4 smaller than the latter (Piran, 2004). However, 
the origin of the variability in the light curves of GRBs remains unclear. 

In this paper, we combine the variability with the screw instability of the magnetic field in the BH 
magnetosphere, and suggest that the variability could be fitted by a series of flares arising from SIBZ, 
which accompanies the release of the energy of the toroidal magnetic field. 

An equivalent circuit MLL'M' for SIBZ is shown in Figure 4a, which consisting of two adjacent 
magnetic surfaces MM' and LL' connecting the BH horizon and the remote load. An inductor is 
introduced in the equivalent circuit by considering that the toroidal magnetic field threads the loop 
MLL'M', and the inductor is represented by the symbol AL in Figure 4a. 

The inductance AL in the circuit is defined by 



where Ll is given by equation (9), and A^/ T is the flux of the toroidal magnetic field threading the 
circuit. The flux A^ T can be integrated over the loop MLL'M' as follows, 



4. 



AN EXPLANATION FOR VARIABILITIES IN GRB LIGHT CURVES 



AL = A^ T /L 



(23) 




(24) 
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where the toroidal magnetic field measured by "zero-angular-momentum observers" is 



B T = 2I L /(aw), (25) 

where a is the lapse function defined in equation (1) (MT82). 

Since the geometric shapes of the magnetic surfaces are unknown, we assume that the surfaces are 
formed by rotating the two radial segments MM' and LL', which span the angle Ad as shown in Figure 
4b. Thus the flux A^f T can be calculated easily by integrating over the region MLL'M'. Incorporating 
equations (23) — (25), we obtain ALas follows. 

I-Lbz [>L bz /M ,~2 , 2 2 a \ 

AL = 2csc0A0/ p 2 dr/A = 2MA6 esc 9 S / \ , % S J df, (26) 

Jr H 7(1+,) (r 2 + a2_2r) 

where f is defined as f = r/M. 

Although the detailed process of SIBZ is still unclear, we suggest that the energy release in one event 
of SIBZ is roughly divided into two stages: the processes for releasing and retrieving magnetic energy, 
respectively. The two processes can be simulated as the corresponding processes in the equivalents — 
Lcircuit. The detailed analysis is given as follows. 

At the first stage the energy of the toroidal magnetic field is released as soon as SIBZ occurs, being 
dissipated on the load and plasma fluid in the way analogous to a discharging process in an equivalent 
R — L circuit. The equation governing the discharging process in R — L circuit is 

dl p 

AL— + {AZ PLSM + AZ L ) I p = 0, (27) 

where AZplsm is the resistance of the plasma fluid in the BH magnetosphere. 

At the second stage the energy of toroidal magnetic field is recovered due to the rotation of the 
BH, and the process for retrieving magnetic energy is modulated by a charging process in an equivalent 
R — L circuit. The equation governing the charging process in R — L circuit is 

dl p 

AL— - + (AZ H + AZ L ) P = Ae H . (28) 
dt 

In equations (27) and (28) AL is the inductance in the circuit MLL'M', and AZplsm is the 
resistance of the plasma in the BH magnetosphere. Incorporating equations (12) and (26), we have 

AL,AZ H = (9.85 x 10- sec) , f""' ^±^M df . (29) 

Combining the initial conditions in the first and second stages, we have the solutions of equations 
(27) and (28) as follows, 
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JP - TP p-t/n (or>\ 
1 disch 1 initial^ ' \ oyJ I 

Ilh = Steady (1 " Z- t/T2 ) ■ (31) 

In equations (30) and (31) If nitial and l v steady are the initial and steady currents, respectively, while 
Idisch an d Ich represent the discharging and charging currents, respectively. The characteristic time 
scales in equations (30) and (31) are given respectively by T\ and r 2 , and they read 

n = AL/(AZ PLSM + AZ L ), (32) 

r 2 = AL/ (AZ H + AZ L ). (33) 
From equations (32) and (33) we have the ratio of T\ to r 2 given by 

n/r 2 = 2AZ h /{AZ PLS m + AZ H ), (34) 

where AZh = AZl is used in deriving equation (34). 

In contrast to disk plasma of perfect conductivity, the resistance AZplsm cannot be neglected 
based on the following considerations: 

(1) The plasma fluid becomes very tenuous after leaving the inner edge of the disk, augmenting 
significantly the resistance due to an increasing radial velocity onto the BH; 

(2) The conductivity of the plasma fluid is highly anisotropic, i.e., the conductivity in the cross-field 
direction is greatly impeded by the presence of the strong magnetic threading the BH (Punsly 2001). 

Although the value of AZplsm is unknown, we can estimate roughly the variability time scales of 
GRBs by combining equation (34) with different cases given as follows. 

CASE I: AZplsm » AZ H leads to t\ « r 2 , and the time scale of two successive flares arising 
from SIBZ is dominated by r 2 . 

CASE II: AZp LS M ~ AZ H leads to T\ ~ r 2 , and the time scale of two successive flares arising 
from SIBZ is about 2r 2 . 

CASE III: AZplsm « AZ H leads to t\ 2r 2 , and the time scale of two successive flares arising 
from SIBZ is about 3r 2 . 

Therefore the variability time scales are insensitive to the values of the unknown AZ PLSM . From 
equation (31) we obtain that the charging current attains 99.3% of l v steady in the relax time t re [ ax = 5t 2 , 
implying the recovery of the toroidal magnetic field, and the variability time scales of GRBs can be 
estimated as follows, 



(St) j = (tsiBz)! ~ 5r 2 , for- AZplsm » AZ H , 



(35) 



(St)jj = (t S iBz) H ~ 10r 2 , for AZ PLSM « AZ H , 



(36) 



(<Jt) m = (t S iBz) ni ~ 15r 2 , /or AZ PL5M « AZ H . (37) 

In CASE I we have T\ « r 2 , implying that the magnetic energy is released much more rapidly 
in the first stage compared with the time for the recovery of the magnetic energy in the second stage. 
Thus CASE I seems more consistent with the feature of the light curves, i.e., an individual pulse is a 
fast-rise exponential decay (FRED) with an average rise-to-decay ratio of 1:3 (Norris et al. 1996). 

By using equations (29), (33) and (35) — (37) we have the variability time scales in the light curves 
of four GRBs in the three different cases as shown in Table 4. In addition, we obtain the curves of (5t) T 
versus a* with the fixed values of n for GRB 990712, GRB 991208 and GRB 021216 as shown in Figure 
5. 

Inspecting Table 4, we find that the variability time scales of tens of msec in the light curves of 
GRBs can be modulated by the two successive flares due to SIBZ, which accompany the BZ process in 
powering the GRBs. We also find that the variability time scales of the four GRBs are generally three 
order of magnitude less than the corresponding durations T 90 , which are consistent with the observations. 

From Figure 5 we find that the curves of (St) T versus a* are almost the same, increasing linearly in 
a rather small slope with the decreasing a*. The values of (St)j remains less than 20 msec during the 
occurrence of SIBZ for these GRBs. 



5. DISCUSSION 

5.1. Mechanism for the recurrent occurrence of SIBZ 

In this paper we discuss the possibility of the modulations of SIBZ on the light curves of GRBs. 
One of the puzzles is what mechanism leads to the recurrent occurrence of SIBZ, and prevents the 
magnetic field from settling to a screw-stable configuration. 

As argued in MT82 the BH magnetosphere consists of a series of magnetic surfaces connecting 
the horizon with the loads, and the total electric current I flowing downward through an m-loop is 
proportional to the toroidal magnetic field B T by Ampere's law. It is argued that these magnetic 
surfaces can be regarded as an equivalent circuit, in which each loop consists of two adjacent magnetic 
surfaces (W02). The total electric current flowing downward through an m-loop is exactly equal to the 
algebraic sum of the poloidal currents flowing in the loops (W03). 

As shown in Figure 1 the critical magnetic surface (henceforth CMS) for SIBZ is represented by 
the critical line MM' . According to the criterion (7) only the toroidal magnetic field outside CMS is 
depressed by SIBZ, while the toroidal and poloidal components of the magnetic field within CMS are 
little affected. On the other hand, the poloidal magnetic field outside CMS still exists in spite of the 
occurrence of SIBZ, and the toroidal magnetic field outside CMS will be recovered because of the twist 
of the poloidal magnetic field arising from the rotation of the BH. So, the rotation of the BH is the 
main mechanism for the recurrent occurrence of SIBZ in the BH magnetosphere. 
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5.2. Rotation period of BH and time scales of recovery of toroidal magnetic fields 

If we take BH mass as 1OM , the rotation period of the BH is only ~lmsec for the BH spin required 
by the criterion of SIBZ. This result implies that toroidal magnetic fields can not be recovered in one 
period of a rotating BH. How to explain the discrepancy between BH rotation and the time scales 
required for recovery of toroidal fields ? 

In spite of lack of detailed knowledge of the screw instability, it is helpful to imagine the magnetic 
field line as an elastic string. The rotating BH always twists the field line, while the field line tries to 
untwist itself. Once the toroidal component of the magnetic field is strong enough to satisfy the criterion, 
the screw instability will occurs, just as a twisted elastic string releases its energy under appropriate 
conditions. In our model we simulate the process for twisting the field line by a transient process for 
accumulating magnetic energy in the inductor AL in the equivalent R — L circuit, which corresponds 
to the increase of the toroidal magnetic field, and the variability time scales of tens of msec in the light 
curves of several GRBs are fitted as the time interval between two successive flares due to SIBZ. Thus 
the threshold of toroidal magnetic field (magnetic energy) cannot be recovered in only one rotation of 
BH, just as a threshold of twisting more than one turn is required for an elastic string to release its 
energy. 



5.3. An explanation for GRBs with XRFs and XRRs 

Recently, much attention has been paid to the issue of X-ray flashes (XRFs), X-ray-rich gamma- 
ray bursts (XRRs) and GRBs, since HETE-2 provided strong evidence that the properties of these 
three kinds of bursts form a continuum, and therefore these objects are probably the same phenomenon 
(Lamb et al. 2004a, 2004b, 2005). The observations from HETE-2 motivate some authors to seek a 
unified model of these bursts. The most competitive unified models of these bursts are off-axis jet model 
(Yamazaki, Ioka, & Nakamura 2002; Lamb et al. 2005) and two-component jet model (Berger et al. 
2003; Huang et al. 2004), in which XRFs, XRRs and GRBs arise from the differences in the viewing 
angles. Unfortunately, a detailed discussion on producing different viewing angles for XRFs, XRRs and 
GRBs has not been given in the above works. 

Our argument on SIBZ and SIMC may be helpful to understanding this issue. It is believed that 
a disk is probably surrounded by a high-temperature corona analogous to the solar corona (Liang & 
Price 1977; Haardt 1991; Zhang et al 2000). Very recently, some authors argued that the coronal 
heating in some stars including the Sun is probably related to dissipation of currents, and very strong 
X-ray emissions arise from variation of magnetic fields (Galsgaard & Parnell 2004; Peter et al. 2004). 
Analogously, if the corona exists above the disk in our model, we expect that it might be heated by the 
induced currents due to SIMC and SIBZ. Therefore a very strong X-ray emission would be produced to 
form XRFs or XRRs. 

Although our model may be too simplified and idealized with respect to the real situation, it 
provides a possible scenario for the occurrence of the screw instability in BH magnetosphere, and it 
may be helpful to understanding some astrophysical observations. We hope to improve our model by 
combining more observations in the future. 
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Fig. 1. — A schematic drawing for magnetic field configuration with screw instability in BH magneto- 
sphere, where the screw-unstable regions for S1BZ and SIMC are indicated by the shaded regions (not 
in scale). 




Fig. 2. — (a) The contours of #s(a*,n) (solid lines); (b) the contours of h c (a*,n) (dashed lines); (c) a 
shaded region for occurrence of SIMC and SIBZ bounded by the contours 9s = and h c = 100. 
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Fig. 3. — BH evolution is described in a* — n parameter space with and without screw instability in 
(a) and (b), respectively. The parameter space is divided by the contours of Pbz = 0, that of P MC = 
and that of A = into several sub-regions. The duration of a GRB is regarded as the time for the 
represent point of BH evolution moving from its initial position in regionIA to the contour Pbz = 0. 
The dimensionless BZ and MC powers are defined as P B z = Pbz/Pq and P M c = Pmc/Pq, respectively. 
The function A = A(a 1r) n) is given in Lei05, being proportional to the rate of change of the BH spin. 
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Black Hole 



Fig. 4. — (a) An equivalent R-L circuit consisting of two magnetic surfaces, MM' and LL', which 
connect the BH horizon and the remote load; (b) a simplified configuration of poloidal magnetic field 
corresponding to the equivalent R-L circuit. 
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Fig. 5.— The curves of (St)j versus a* with n =3.54, 3.55 and 3.58 for GRBs 990712, 991208 and 991216 
in (a), (b) and (c), respectively. 
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Table 1. The Main Features of Several Models Invoking the BZ Process For Powering GRBs 



Model 


BZ Process 


MC Process 


Surrounded By 


Objects 


Half-Opening Angle 


Variability 


I 


Yes 


No 


No disk 


GRB 


No 


No 


II 


Yes 


No 


Transient Disk 


GRB 


No 


No 


III 


Yes 


No 


Torus 


GRB 


No 


No 


IV 


Yes 


Yes 


Disk 


GRB/SN 


No 


No 


V 


Yes 


Yes 


Torus 


GRB/SN 


Yes 


No 


VI 


Yes 


Yes 


Disk 


GRB/SN 


Yes 


No 


This Model 


Yes 


Yes 


Disk 


GRB/SN 


Yes 


Yes 



Table 2. Five GRBs of True Energy E 1 and T 90 s Fitted with Different Power-law Index n and B H 

with the Predicted E S n- 



GRB a 




T90 


n 




p>h 








970508 


0.234 


15 d 


3.89 


3.46 


0.97 


0.70 


1.947 


1.732 


990712 


0.445 


30 e 


3.98 


3.54 


0.81 


0.61 


1.989 


1.816 


991208 


0.455 


39.84' 


3.99 


3.55 


0.75 


0.56 


1.995 


1.817 


991216 


0.695 


7.51' 


4.06 


3.58 


1.80 


1.40 


2.032 


1.887 


021211 


0.117 


13.3 s 


3.81 


3.42 


0.9 


0.6 


1.901 


1.681 



Note. - 

a This sample of GRB-SNc is taken from Table 1 of Dar (2004). 

''The true energies E~ of the GRBs are given in terms of 10 51 ergs from 

Table 1 of Frail et al (2001). 

c The predicted energies Egpf of the SNe are given in terms of 10 51 ergs 
based our models. 

d The duration of GRB 970508 is from Costa ct al. (1997). 

e The duration of GRB 990712 is from Heise et al. (1999). 

'The durations of GRBs 991208 and 991216 are from Table 2 of Lee & Kim 

(2002). 

9 The duration of GRB 021211 is from Sakamoto et al. (2005). 
h Bi§ is the magnetic field at the BH horizon in terms of 10 15 gauss. The 
initial BH spin is assumed to be a*(0) = 0.9, and the data in the left and 
right sub-columns of the parameters n and Bis and Egpf are calculated 
based on Model VI and this model, respectively. 



Table 3. Five GRBs Fitted with Different fsz and fuc Invoking Model VI and This model. 



GRB 


a? RB 




Sbz 




Sue 








970508 


0.798 


0.871 


0.009 


0.020 


0.991 


0.980 


0.233 


0.226 


990712 


0.767 


0.853 


0.014 


0.028 


0.986 


0.972 


0.234 


0.227 


991208 


0.763 


0.852 


0.014 


0.029 


0.986 


0.971 


0.235 


0.227 


991216 


0.736 


0.842 


0.019 


0.032 


0.981 


0.968 


0.236 


0.228 


021211 


0.823 


0.880 


0.006 


0.016 


0.994 


0.984 


0.232 


0.225 
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Table 4. Variability Time Scale in the Light-Curves of Five GRBs. 



GRBs 




Parameters 






Tgo(s) 




Variablity Time scale (ms) 






a* 


n 


r s 


h c 




(St) j 


(St)ii 


Win 


970508 


0.9 


3.46 


5.96 


87.71 


15 


13.42 


26.84 


40.25 


990712 


0.9 


3.54 


6.45 


85.56 


30 


13.23 


26.46 


39.68 


991208 


0.9 


3.55 


6.52 


87.71 


39.84 


13.42 


26.84 


40.25 


991216 


0.9 


3.58 


6.77 


87.70 


7.51 


13.42 


26.83 


40.25 


021211 


0.9 


3.42 


5.77 


59.32 


13.3 


10.90 


21.80 


32.70 



Note. - 

rs = r s/M is the critical dimcnsionlcss radius for SIMC. 



